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ABSTRACT 



Lyman break galaxies (LBGs) represent one of the kinds of star-forming galaxies that are found in the high-redshift universe. The 
detection of LBGs in the FIR domain can provide very important clues on their dust attenuation and total SFR, allowing a more detailed 
study than those performed so far. In this work we explore the FIR emission of a sample of 16 LBGs at z ~ 3 in the GOODS-North and 
GOODS-South fields that are individually detected in PACS-lOOyum or PACS-160yum. These detections demonstrate the possibility 
of measuring the dust emission of LBGs at high redshift. We find that PACS-detected LBGs at z ~ 3 are highly obscured galaxies 
which belong to the Ultra luminous IR galaxies or Hyper luminous IR galaxies class. Their total SFR cannot be recovered with the 
dust at tenuation factors obtained from their UV continuum slope or their SED-derived dust attenuation employing 'Br uzual & CharJod 
( 120031) templates. Both methods underestimate the results for most of the galaxies. Comparing with a sample of PACS-detected LBGs 
at z ~ 1 we find evidences that the FIR emission of LBGs might have changed with redshift in the sense that the dustiest LBGs found 
at z ~ 3 have more prominent FIR emission, are dustier for a given UV slope, and have higher SFR for a given stellar mass than the 
dustiest LBGs found at z ~ 1. 

Key words, cosmology: observations - galaxies: stellar populations, morphology, infrared. 



1. Introduction 



X 



The dropout technique, that segregates the so-called Lyman 
break galaxies (LBGs), is one of the most employed and suc- 
cess ful methods to find high-red s hift star- forming ( SF) galax- 
ies ("Burgarella et a l.ll2()06tl20q7t IChen et"aLll2013l: lOteo et al.' 
2013; Basu-Zvch e t all 120111: iHathi et all l2012t iM adau et al. 
1996; Steideletal. 1996, 1999, 2003; Stanwavetal. 2003; 
Giavalisco et al.. 2004: .Bunker et al.. .2004; .Verma et al.. .2007; 
Iwata et aP l2007h . Many studies have analyzed the physical 
properties of LBGs with photometric and spectroscopic data, 
from X-ray to radio wavelengths. The detection rate of LBGs 
is very dependent upon their redshift and the wavelength of 
the observations. LBGs can be easily detected in the opti- 
cal with the current deep cosmological surveys. A fraction of 
them are detected near-IR emission. The situation becomes 
problematic when we move to redder wavelengths and study 
high-redshift (z > 3) LBGs. Only massive LBGs at z ~ 3, 
log( M/Mo) ~ 11, are detected in IRAC or MIPS channels 
(Mag dis et ai]l2010ah . The detection of some LBGs at diff'er- 
ent redshifts with MIPS-24yum revealed a population of lumi- 
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* Herschel is an ESA space observatoiy with science instruments pro- 
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nous infrared LBGs (Huang et al. 2005; RigoDoulou et al.'2006t 
iBurg arella et al. 2006. 2007: Basu-Zvch et al. 201 1). On the FIR 
side,[Burgarella et all (1201 ll) found SPIRE detections for a sam- 
ple of 12 GALEX-selected LB Gs at z ~ 1 and one LBG at 
z ~ 2. iRigopoulou et al.l (1201 Ol) did not found individually de- 
tected LBGs at z ~ 3 in SPIRE bands and lMagdis et alj <2010bl) 
did not found PACS-FIR counterparts of their LBGs at z ~ 3 
either. In redder wav elengths than FIR, very few LBGs have 
been directly detected ("Chapman et al."200C| IChapman & CasevI 
l2009l) . Recendy, Magdis et al. (2012) present CO[3^2] obser- 
vations of two PACS-detected FIR-bright LBGs at z ~ 3 which 
indicate that the steep evolution of Mgas/M» of normal galaxies 
up to z ~ 2 is followed by a flattening at higher redshifts, pro- 
viding evidence for the existence of a plateau in the evolution 
of the specific star f ormation rate (sSFR) at z > 2.5 (see also 
lBouweiiseTani2012h . 

In order to understand high-redshift LBGs a deeper study 
of their FIR spectral energy distribution (SED) is needed. In this 
work we exploit the deep publicly available FI R data taken in the 
framework of the GOODS-Herschel project (lEIbaz et al .11201 ih 
in the GOODS-North and GOODS-South fields to increase the 
sample of FIR-bright PACS/SPIRE-detected LBGs at z ~ 3 and 
analyze in more detail their FIR emission. Additionally, we com- 
pare the behavior of FIR-bright LBGs at z ~ 3 with their ana- 
logues at z ~ 1. This letter is structured as follows: the selection 
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Fig. 1. Differential SED-derived properties between PACS-detected (sliaded green histograms) and PACS-undetected (grey histograms) LBGs 
at z ~3. T he SED-derived properties have been derived by fitting the observed rest-frame UV to near-IR photometry of the galaxies with 
iBruzual & Chariot (2003) templates associated to a constant SFR and fixed metallicity of Z = 0.4Zq. Histograms have been normalized to 
their maxima in order to clarify the representations. 



of our LBGs and their FIR emission is explained in Section |2] 
In Sections[3]and|4]we present and discuss the properties of our 
PACS-detected galaxies. Finally, the main conclusions of the 
work are summarized in Section |5] Throughout this letter we 
assume a flat universe with (Q,„, Oa, ho) - (0.3, 0.7, 0.7) an d all 
magnitudes are listed in the AB system (lOke & Gunnlll983h . 

2. FIR emission of LBGs at z ~ 3 

We select our LBGs at z ~ 3 by employing the classical dropout 
technique with the broadband filters U, V, and /. Details of 
the analytical selecti on criterion can be found in Appendix |A] 
Briefly, a large set of (iBruzual & Chariotll2003l hereafter BC03) 
templates at different redshifts were convolved with the trans- 
mission curves of the U, V, and / filters. Analytically, we select 
our LBGs through: 



U -V > l.llxiV -I) + Q.51 



(1) 



With this, we isolate a sample of 2652 and 2537 LBGs can- 
didates in GOODS-North and GOODS-South, respectively. We 
look for their FIR counterparts by employing PACS and SPIRE 
publicly available data coming from the GOOT)Si-Herschel 
project. In this step we adopt a matching radius of 2". We ob- 
tain that, among the whole sample of 2652 and 2537 LBGs in 
GOODS-North and GOODS-South, 7 and 9, respectively, are 
likely detected in any of the PACS bands. These numbers repre- 
sent percentages of FIR detection lower than 0.5%. One LBG in 
GOODS-South and another in GOODS-Noith are spectroscop- 
ically confirmed to be at z ~ 3 (see Figure IaTTT i. Additionally, 
among the 7 PACS-detected LBGs in GOODS-North, 5 of them 
are also detected in SPIRE-250yum, SPIRE-350yum, or SPIRE- 
500yum. No SPIRE data is publicly available in GOODS-South 
by the time of writing of this paper and, therefore, we cannot 
report on SPIRE-detections in that field. By using 15"xl5" 
optical and mid-IR cutouts of our PACS-detected sources we 
have checked that source confusion is unlikely in our sample 
(see Figure lATJb . Details on rest-frame UV-to-near-IR and FIR 
SED fitting can be found in Appendix |A] We obtain the total 
SFR of our PACS-detected LBGs by combining their UV and 
FIR measurements, S FRtoVdi = •S/^/^uv.unconected +S FRm, where 
SFRuYiinconected IS the SFR associatcd to the dust-uncoiTected 
rest-frame UV luminosity and S FRir is the SFR associated to 
the IR emissi on of our galaxie s. Both terms are calculated by 
employing the 'Kennicut^ ( Il998l) calibrations. 

Figure [1] indicates that PACS-detected LBGs are more mas- 
sive (median log (M./Mq) = 10.7 + 0.8), have higher SED- 
derived dust attenuation (median E^iB - V) = 0.40 + 0.15), 



stronger Balmer 4000A break (median Lr/Ll = 1-44 + 0.21), 
and redder UV continuum slope (median j3 - -0.74 + 0.56) than 
PACS-undetected LBGs. The uncertainties indicate the width 
of the distributions. There is no significant difference between 
the age of both populations. Table IaTI compiles some proper- 
ties of our PACS-detected LBGs at z ~ 3. The median red- 
shift of the sources is Zphot = 3.1. It can be seen in Table IaTI 
that all our PACS-detected LBGs are Ultra luminous IR galax- 
ies (ULIRGs) or Hyper luminous IR galaxies (HyLIRGs). This 
is a direct consequence of the depth of the FIR data employed: 
only the FIR-biightest sources are detected. The IR/UV-derived 
dust attenuation is larger than 5 mag in 1200 A for most of the 
galaxies. This indicates that, despite being selected through rest- 
frame UV color and being UV-bright galaxies, it exists subpopu- 
lation of LBGs at z ~ 3 that is very dusty. PACS-detected LBGs 
are extreme SF galaxies with UV-i-IR-derived total SFR typically 
higher than 500 Mq yr"' . 

We compare in the next sections the FIR-derived properties 
of our PACS-detected LBGs at z ~ 3 with those of a sample 
of 32 PACS-detected LBGs at z ~ 1 in the GOODS-North and 
GOODS-South field (see Appendix lAl for details on the sample 
selection). There is a fundamental difference between the two 
populations of FIR-bright LBGs: at z ~ 1, all PACS-detected 
LBGs are LIRGs with log(LiR/L0) < 11.7 whereas all PACS- 
detected LBGs at z ~ 3 are ULIRGs with log(LiR/Lo) > 12.4 
or HyLIRGs. LIRGs are not detected at z ~ 3 due to the FIR 
selection bias but ULIRGs at z ~ 1 could have been found since 
the observations are complete down to the expected brightness 
for ULIRGs in this redshift range. The fact that we do not find 
any ULIRG-LBGs with log(LiR/L0) > 12.4 at z ~ 1 might in- 
dicate that the FIR emission of LBGs might have changed with 
redshift. When analyzing the evolution of the FIR emission of 
LBGs with redshift, the difference in the surveyed volumes need 
to be taken into account. Both LBGs at z ~ 1 and at z ~ 3 are 
located in the same fields and the comoving volume within the 
redshift range 2.5 < z < 4 is about 5.4 times higher than that in 
0.8 < z < 1.2. Consequently, the fact that we do not find any 
ULIRG with log (Lir/Lq) > 12.4 in the sample of LBGs at z ~ 1 
might be also due to the smaller volume surveyed. However, in 
Oteo et al. (2013b, submitted) we do not find any LBG with an 
ULIRG nature at z ~ 1 either in a volume which is only two 
times smaller than that surveyed in the present work at z ~ 3. 
Supposing that the density of ULIRGs with log (Lir/L©) > 12.4 
is the same at z ~ 3 and z ~ 1, we should have found 9 galaxies 
of that type at z ~ 1 in Oteo et al. (2013b). The fact that all 
of the PACS-detected LBGs at z ~ 1 have log(LiR/L0) < 11.7 
reinforces the idea of an evolution of the FIR emission of LGBs 
with redshift. This would be in agreement with the evolution of 
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Fig. 2. FIR-derived properties of our PACS-detected LBGs at z ~ 3 (light green dots for galaxies in GOODS-North and dark green dots for 
galaxies in GOODS-South) and at z ~ 1 (orange open squares). The grey shaded zones represent the accessible areas with PACS at z ~ 3 for 
log (Luv/^g) > 10-0 or for log (Luv/^o) > 1 1.0, which represent the range of the rest-frame UV luminosity for LBGs at z ~ 3. In the bottom left 
and bottom right panels, blue open triangles represent the whole sample of LBGs at z ~ 3 in GOODS-South regardless the detection in PACS. The 
patterns seen in the blue triangles are due to the sampling of the BC03 templates employed to fit the rest-frame UV-to-near-IR SED. In the fourth 
panel, the red straight line is the linear fit to the whole population of LBGs at z ~ 3 and the red dashed line is a linear fit to the PACS-detec ted 
LBGs at z ~ 3. Upper left: Dust attenuation as parametrize d by the rat i o betw een the total IR a nd rest-frame UV lum inosities jBuat et al.l2005l) as 
a function of the UV continuum slope. The relations of .Meurer et alj ( 11999 ) for local SB andlBoissier et al.l (l2007h f or loca l normal SF galaxies 
are represented with dashed and continuous curves, respectively. We also plot the coiTection of Takeuchi e t al.l (120121) to the lMeurer et al.l ( 119991) 
relation with a dashed-dotted curve. Upper right: UV-l-IR-derived total SFR against SED-derived dust corrected total SFR. Solid line is the one- 
to-one relation. Bottom left: SFR versus stellar mass plane. The solid line is the MS of galaxies at z ~ 2 reported in Daddi et al. (2007) . Th e 
dashed and dashed-dotted lines are 4 and 10 times the MS and are employed to separate normal SF galaxies from SB (iRodighiero et al.ll201lh . 
Bottom right: sSFR versus stelar mass plane. 



the evolution of the SFR density (SFRD) of the univers e between 
z ~ 1 andz ~ 3 (see for example lHopkins & Beacoml2006.) . The 
higher values of the SFRD at z ~ 3 might favor the appearance 
of such as extreme LBGs at high redshift. 



3. Dust-correction factors 

The upper left panel of Figure ^indicates that, for each value of 
the UV continuum slope (J3), the dust attenuation parametrized 
by IRX = log (Lir/Luv) of the PACS-detected LBGs at z ~ 
3 is much higher than that predicted by the star-burst (SB) 



relations of iMeurer et all (Il999h and iTakeuchi etaH (l2012l) . 
This is in agreement with FIR-bright galaxies at low redshift 
(iGoldader et al.ll2002l) . Therefore, there is a population of LBGs 
at z ~ 3 for which the SB relations fails at recovering their dust 
attenuation from their UV continuum slopes. Specifically, this 
method would produce underestimated results. For a given UV 
continuum slope, the dustiest LBGs at z ~ 1 have lower dust at- 
tenuation than the dustiest LBGs at z ~ 3. For each /?, LBGs at 
z ~ 1 as dusty as those at z ~ 3, could have been found since the 
observations are complete. The absence of this kind of sources 
suggests that the upper envelope of the locus of LBGs in a IRX-yS 
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diagram might have changed with redshift. Interestingly, despite 
PACS -detected LBGs at z ~ 3 are dustier than those at z ~ 1, the 
UV continuum slopes of both populations span within a similar 
range. 

The upper right panel of Figure |2] indicates that the SED- 
derived dust attenuation of our PACS-detected LBGs at z ~ 3 
cannot be used to recover their U V +IR-derived total SFR. The 
SED-derived total SFR is underestimated for most of the galax- 
ies. This might be a consequence of a clumpy or patchy ge- 
ometry of the dust regions, and the UV-luminous part and IR- 
luminous part are different in each galaxy. The reported under- 
estimation is in agreement with previous results suggesting that 
SED-fitting procedures fails for e stimating the total SFR in dust- 
obscu red high-redshift galaxies (lOteo et al.ll2012t IWuvts et all 
1201 lb . At z ~ 1, instead of an underestimation, the SED-derived 
dust attenuation overestimates the total SFR for LBGs. The 
SED-derived total SFR of the PACS-detected LBGs at z ~ 1 and 
z ~ 3 spans over a similar range, whereas the UV-i-IR-derived 
SFR is much higher for PACS-detected LBGs at higher redshifts. 
The total SFR might be also recovered from the UV continuum 
slope by applying, for e xample, th e Meureret al. (1999) law or 
the correction given by [T akeuchi et al. (2012) ( Oteo et al.ll2013l: 
iBouwens et alj|2009ll201 L.Castellano et al...2012l) . In this case, 
since our PACS-detected LBGs are well above the SB relations, 
this technique would also underestimate the real dust attenua- 
tion, and thus the total SFR, of the studied galaxies. Therefore, 
we conclude that the only way to obtain reliable values of the 
dust attenuation and total SFR of our PACS-detected LBGs at 
z ~ 3 is using their UV and FIR emission. This highlights the 
importance of FIR data when studying high-redshift sources. 

4. SFR and stellar mass 

The location of our PACS-detected LBGs at z ~ 3 in an SFR 
versus stellar mass diagram is shown in the bottom left panel of 
Figure |2] For comparison, we also plot the location of the whole 
population of LBGs at z ~ 3 in GOODS-North regardless their 
detection in the FIR. For a given stellar mass, PACS-detected 
LBGs have higher total SFR than those PACS-undetected. It 
should be remarked that the total SFR for PACS-detected and 
PACS-undetected galaxies has been calculated differently. For 
the former, the total SFR is derived from their direct UV and FIR 
emission. For the latter, the total SFR has been obtained by cor- 
recting the rest-frame UV luminosity with the best-fitted value 
of the SED-derived dust attenuation, as explained in Section |2] 
We also show in that panel the SFR-stellar mas s relatio n of the 
main sequ ence (MS) at z ~ 2 rep orted in Dad di et aP (|2007). 
Following Rodighiero et al. (201 1) we represent 4 and 10 times 
the MS of Daddi et al. (2007 ) as a r eferen ce for the definition 
of SB galaxies (see also Elb az et all (1201 ih ). At z ~ 3, the low 
number of FIR detections of SF galaxies has prevented the defi- 
nition of an MS. Thus we have to base our results on the MS at 
z ~ 2 without forgetting that there are evide nces of an evolution 
of the MS with redshift (lElbazet al.ll2011h . It can be seen that 
most PACS-detected LBGs at z ~ 3 are located above 4 times the 
MS and, therefore, have a SB nature. The FIR-brightest LBGs 
at z ~ 3 have higher values of the total SFR for a given stellar 
mass than the FIR-brightest LBGs at z ~ 1 . This indicates an 
evolution of the high-SFR tail of LBGs with redshift. 

The sSFR versus stellar mass for our PACS-detected LBGs 
at z ~ 3 is presented in the bottom right panel of Figure |2] It 
can be seen that PACS-detected LBGs at z ~ 3 follow a linear 
relation as a direct consequence of the limiting FIR fluxes of the 
observations employed: only very FIR-bright galaxies with high 



total SFRs are detected. For each stellar mass, PACS-detected 
LBGs at z ~ 3 have higher values of the sSFR than those PACS- 
undetected. For a given stellar mass, the FIR-brightest LBGs at 
z ~ 1 have lower values of the sSFR than those at z ~ 3. This 
again indicates an evolution of the FIR emission of LBGs with 
redshift, at least in the FIR-brightest side. 

5. Conclusions 

In this work we have found a sample of 16 FIR-bright LBGs 
at 2.5 < z < 4.0 in GOODS-North and GOODS-South fields. 
These galaxies are individually detected in PACS-100/im or 
PACS-160yL(m, probing directly their dust emission. These detec- 
tions allow us to determine their total IR luminosities, dust atten- 
uation, and total SFR without the uncertainties that introduce the 
traditionally used SED-fitting technique with BC03 templates. 
The dust attenuation and total SFR of these objects cannot be 
recovered from the dust correction factors obtained with their 
UV slope or their SED-derived dust attenuation since both meth- 
ods underestimate the results. The only way of obtaining ac- 
curate results is employing UV and FIR data, highlighting the 
importance of these wavelengths for the studies of high-redshift 
sources. Comparing with a sample of PACS-detected LBGs at 
z ~ 1 we find evidence that the FIR emission of LBGs might 
have evolved with redshift: the dustiest LBGs at z ~ 3 have 
more prominent FIR emission, have higher dust attenuation for 
a given UV slope, and have higher total SFR for a given stellar 
mass than the dustiest LBGs at z ~ 1 . 
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Fig. A.l. Left plot: location of the LBGs of ISteidel et alJ i2003h with confirmed spectroscopic redshifts (red dots) in the color-color diagram 
employed to select the LBGs studied in this work. The light grey sh aded zone represen ts our selection window for LBGs at z ~ 3 (see Equation 
lA.lt . Black open squares are th e complete sample of galaxies in the lCapak et al.l l l2004h photometric catalog. For comparison, we also show the 



selection window for LBGs of Pentericci 

z 



e sarn p 
eTaLl d: 



120101) with a dark grey shaded zone. Right plot: Photometric redshift accuracy for our LBGs at 
3. Only galax ies with confirmed spectroscopic redshift via emission lines (light green symbols) or absorpti on lines (dark green symbols) in 
ISteidel et alj l l2003l) are included. The PACS-detected LBG spectroscopically confirmed to b e at z = 2.92 9 in Steidel et al. ( z003) is indicated 
with a orange dot. The PACS-detected LBG spectroscopically confirmed to be at z = 2.975 in lPopesso et al, (,2009,) and ,Balestraet al., (,2010) is 
presented with a red dot. 



Appendix A: Sample selection and SED fitting 

Our LBGs are selected through t he classical dropout t echniq ue with the U, V, and / filters. By convolving the transmission curves 
of those filters with a large set of iBruz ual & Chariot (2003) (hereafter BC03) templates associated to different physical properties 
and redshifts (0 < z < 6) (see for example .Oteo et al.,,2013.) . we derive the following analytical selection criterion; 



t/-y> 1.71 x(y-/)-H 0.57 (A.l) 

Th e mu lti-wavelength ph otometric information from the optical to mid-IR is taken from 'Sant ini et al.l (l2009h . ICapak et aP 
(|2004 . and lWang et al.l (l2()10i) . We focus our work on the GOODS-North and GOOD S-So uth fields since these are the two fields 
covered with the GOODS-Herschel FIR obser vations. The criterion shown in EquationlA.ll is less restrictive than others employed 
in previous studies to look for LBGs at z ~ 3 (iGrazian et a l."2007': Pentericci et al. 2010). It can be seen in the left panel of Figure 
A. II that our selection window covers a wider area in the color-color diagram than that represented by the JPentericci et al. (120 lOi) 
criterion. We also show in the left panel of Figure lA. li the location of the spectroscopically confirmed LBGs at z ~ 3 of lSteidel et al.l 
(l2003l) in the GOODS-North field. T o do that w e have used the same photometric information than that employed to look for our 
LBGs, i.e. the photometric catalog of ICapak et a l. (2004). Most of the spectroscopically confirmed LBGs at z ~ 3 of lSteidel et al.l 
(2003) are within our selection window and, therefore, we conclude that Equation I A. 1 1 truly segregate LBGs at z ~ 3. In order to 
clean the sample from lower-redshift interlopers, we only select those galaxies whose photometric redshifts (see next paragraph) are 
within 2.5 < Zphot ^ 4. We discard AGNs by ruling out sources with X-ray emission. In total, we segregate 2652 and 2537 LBGs 
candidates in GOODS-Noith and GOODS-South, respectively. 

Photometric redshifts, age, dust attenuation, and stellar mass of our LBGs are obtained simultaneously by fitting their U to 
IRAC-8jum with a large set of BC03 templates associated to different physical properties of galaxies. LBGs at z ~ 3 are not detected 
in the UV and the MIPS-24yum have a significant contribution of dust emission features which are not considered in the elaboration 
of the BC03 templates. For these reasons, these wavele ngths are not em ployed in the SED fits. We build the BC03 templates by 
using the software GALAXEV. In this process we adopt a'Salpeterl (1 19551) initial mass function (IMF) distributing stars from 0.1 to 
100 M0 and select a constant value for metallicity of 2=0.42©. We consider values of age from 1 Myr to 7 Gyr in st eps of 10 Myr 
from 1 Myr to 1 Gyr and in steps of 100 Myr from 1 Gyr to 7 Gyr Dust attenuation is included in the templates via the lCalzetti et al.l 
(I2OOOI) law and parametrized through the color excess in the stellar continuum, Es(B - V). We selec t values for E/B - V) ranging 
from to 0.7 in steps of 0.05. We include intergalactic medium absorption adopting th dMadaul ( ll995h prescription. Regarding SFR, 
we adopt time-constant models. In this case, different values of the SFR does not change the shape of the templates, and the SFR 
can be obtained by using the Kennicutt ( 1998i) calibration: 



SFRr 



,rrecUMoyr-'] = lAxlO-"'L 



(A.2) 



where L1500 is the rest-frame UV luminosity in 1500A. The Lisoo is obtained for each galaxy by convolving its best-fitted template 
with a top-hat filter centered in rest-frame 1500 A. It should be noted that, throughout the work, we distinguish between Lj/y 
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defined in a vL,, way and L1500 considered in Ly units. The SFR derived from Equation lA. 2 1 is uncorrected for the attenuation that 
dust produces in the SED of galaxies. In order to obtain an estimation of the dust-corrected total SFR we have to introduce into 
Equation lA.2l the dust-corrected Lisoo- It is obtained from L1500 by multiplying it by the dust correction factor lO""^'*'™", where A1500 
is the dust attenuation in 1500A. The values of A 1500 are obtained from the SED-derived Es(B - V) assuming the ICalzetti et al.l 
(I2OOO ) law. Throughout the work, the total SFR calculated in this way will be called SED-derived total SFR, in contraposition with 
the more accurate UV+IR-derived total SFR obtained with the direct emissions in the UV and FIR (see later in the text). Once both 
age and dust-corrected total SFR are known for each source, and according to the assumed time-independent SFH, the stellar mass 
can be obtained from the product of both quantities. 

We define the amplitude of the Balmer 4000 A break as the ratio between the rest-frame 4500 and 3500 A luminosities. These 
are also calculated by convolving each best-fitted template with top-hat filters centered in 4500 and 3500 A, respectively. The UV 
continuum slopes are calculated by fitting the UV continuum of each best-fitted template with a power law function (see for example 
lOteo et al]|2013l: iFinkelstein et al.1 20111) between rest-frame 1300 and 3000 A. We have not included the contribution of emission 



Hnes in the SEP fits (IZackrisson et al.ll2008t ISchaerer & de Barrosll2009l 120101: IS chaerer et al."201 ll 1201 3l Ide Barros et aDl2012h 



since all our R4CS-detected LBGs have mae^m - ?M4 5^in > and, according to Ide B arros et al. (2012), the contribution of emission 



lines in their SED-derived parameters is not expected to be significant. The right panel of Figure lATI s hows the a ccuracy of the 
photometric redshifts of our LBGs, that we define as cr az - \Zphot - Zspec 1/(1+ Zspec) (see for example,Oteo et al.ll2013i) . is lower than 



0.2 for all the sources. These values are comparable with those obtained at lower redshifts (IHaberzettl et al.ll2012l) and it is enough 
for our purposes. 

For all our PACS-detected LBGs we fit ICharv & Elba j (1200 ll) templates to their IRAC-8/im, MIPS-24/.<m, PACS, and SPIRE 
fluxes (when available). As an example of the typical FIR SED-fitting results, we show in Figure IA!2] the rest-frame UV-to-FIR 
SEDs of the 7 PACS-detected LBGs at z ~ 3 in GOODS-North. We then obtain their total IR luminosities by integrating each 
best-fitted template in the rest-frame interval [8-1000] ;um. The dust attenuation of the PACS-dete cted LBGs is parametrized by 
the ratio between the total IR and rest-frame UV luminosities (see foe example iBuat et al]|2005l) . Their total SFR, SFRxaiaX - 
SFR\jY + SFRiR , are calculated by combining the total IR and rest-frame UV luminosities and applying the Kennicutt (1998) 
calibrations: SFRuv,uncoiTected[A^0yr"'] = 1-4 x lO^^'^Lisoo and SFRiR[M0yr"'] - 4.5 x 10"'*'*Lir. It should be remarked that this 
procedure for obtaining the total SFR assumes that all the light absorbed by dust in the rest-frame UV is reemitted in turn in the FIR 
(see for example Magdi set al.ll2010ah . 

For comparison, and with the aim of carrying out evolutionary studies, we compare the FIR-derived properties of our PACS- 
detected LBGs at z ~ 3 with a sample PACS-detected LBGs at z ~ 1 in the GOODS-North and GOODS-South fields. At that 
redshift, the Lyman break is located in the UV and LBGs must be se lected by employi ng UV colors from space-based observations. 
For segregating our LBGs at z ~ 1 we adopt the selection criterion of lOteo et al.l(l2013h . The photometric redshifts, stellar mas s, age, 
and dust attenuation of the galaxies are obtained from a SED-fitting procedure with BC03 to their UV to near-IR fluxes (Capa ket al.l 
|2007) in the same way that it was done for LBGs at z ~ 3. Their FIR emission is also characterized with PACS data taken from 
the GOODS-He rschel project. The FIR-derived properties of the PACS-detected LBGs at z ~ 1, i.e. dust attenuation and UV-i-IR- 
derived total SFR are obtained in the same way than it was done for the PACS-detected LBGs at z ~ 3. Due to the observational 
bias, LBGs at higher redshifts tend to be brighter Therefore, if we want to compare galaxies at different redshifts which are selected 
through the same selection criterion, we must limit both samples to the same rest-frame UV luminosity. At z ~ 3, most LBGs have 
log Luv/Lq > 10. Imposing this limit, we end up with a sample of 3 1 PACS-detected LBGs at z ~ 1 . 



Table A.l. Proper ties of the PACS- detected LBGs at z ~ 3 studied in the pr esent work. We show the galaxy ID, sky position, photometric redshift, total IR luminosity, dust attenuation derived from 
the Lir/Luv ratio JBuat et al.l2005h . UV+IR-derived SFR ( lKennicuttlll998l) . stellar mass, and UV continuum slope. 



ID 



RA[deg] Dec[deg] Zphot log(LiR/LB) Ai2oo[mag] ^ FJ^totai [Mp yr ' ] log (MJMq) UV slope 
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Fig. A.2. Rest-frame UV-to-FIR SEDs of the 7 PACS-detected LBGs at z ~ 3 in the GOODS-North fields. These plots represent the typical 
SED of the studied galaxies. O range curves represent the best-fitted Bmzual & Chariot (2003) templates to the U-band to IRAC-8yum fluxes. Red 
curves represent the best-fitted ICharv & Elbazl(l200lh templates to the MIPS-24yum to PACS and SPIRE fluxes. The photometric redshift of each 
source is indicated in each panel. 
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Fig. A.3. ACS optical cutouts (15"xl5") of the nine PACS-detected LBGs located in the GOODS-South field. In each image, the PACS- 
detected LBG is the source located in the center. Blue contours represent the emission in the bluest available IRAC image. Red contours trace the 
MIPS-24yum emission. We use MIPS-24yjm instead of PACS countours because they are visually clearer and the PACS fluxes are extracted with 
MIPS-24y[im priors. 



